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Abstract 

We review tiie current issues of nonequilibrium phase transitions, in particular, in the early uni- 
verse. Phase transitions cannot maintain thermal equilibrium and become nonequilibrium when 
the thermal relaxation time scale is greater than the dynamical time scale. Nonequilibrium phase 
transitions would have happened in certain evolution stages of the early universe because the 
rapid expansion quenched matter fields. We apply the recently introduced Liouville-von Neumann 
method, another canonical method, to nonequilibrium phase transitions in the Minkowski space- 
time and find the scaling behavior of domain sizes. Topological defects are thus determined by the 
dynamical processes of nonequilibrium phase transitions. However, the expansion of the universe 
freezes domain growth in the comoving frame and thus leads to a scale-invariant domain size. We 
also discuss the physical implications of nonequilibrium phase transitions in the early universe. 
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I. INTRODUCTION 



A system undergoes a phase transition when its symmetry is broken exphcitly [1, 2, 3]. 
The phase transition proceeds either in equihbrium or in nonequihbrium (out of equihbrium) 
depending on the ratio of the thermal relaxation time to the dynamical time. When the 
thermal relaxation time scale is shorter or longer than the dynamical time scale, the system 
evolves in equilibrium or out of equilibrium. In particular, a quenched system undergoes 
nonequihbrium phase transition when the quench rate is faster than the relaxation rate. 
Matter fields are believed to have undergone such nonequihbrium phase transitions in the 
early universe as the universe expanded and thereby temperature dropped rapidly. It is 
likely that such nonequihbrium phase transitions of matter fields can also be realized and 
observed in RHIC and LHC experiments in the near future. 

The finite-temperature field theory has been the most popular approach to equilibrium 
phase transitions [4]. The effective potential of quantum fluctuations around a classical 
background provides a convenient tool to describe phase transitions. However, quantum 
fluctuations of long wavelength modes suffer from instability during the phase transition, 
become unstable and grow exponentially. This is the origin of the imaginary part of the 
effective action for phase transitions, which gives the decay rate of the false vacuum [5]. Thus 
the finite-temperature effective action does not properly take into account the dynamical 
processes of phase transitions. 

On the other hand, nonequihbrium phase transitions have been frequently treated in the 
closed-time path integral defined in a complex time plane by Schwinger [6] and Keldysh 
[7]. Other methods are the Hartree-Fock or mean field method [8, 9, 10, 11, 12], the 
1/A^-expansion method [13, 14, 15], the time-dependent variational principle [16], and etc. 
In this paper we use the recently developed Liouville-von Neumann (LvN) method [17, 
18, 19, 20, 21]. The LvN method is a canonical method trying to solve the functional 
Schrodinger equation [22, 23]. The quantum LvN equation provides invariant operators 
that solve exactly the functional Schrodinger equation [24] . Compared with other methods, 
this canonical method has advantages that it can be applied directly to time-dependent 
systems and density operators are readily found. The mode-decomposed Hamiltonian of 
quenched models for phase transitions has this time-dependent feature. 

We use quenched models of theory with the sign changing mass term for nonequilib- 
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rium phase transitions. These quenched models have quench rates that determine exphc- 
itly time-dependent phase transition processes, an instantaneous quench model [9, 10, 18] 
and a finite quench model [18, 25]. In both models the field, after symmetry breaking 
phase transition, begins to roll from the false vacuum toward the true vacuum. During 
this period of phase transition, the system undergoes spinodal instability and its long 
wavelength modes grow exponentially. Domains grow through spinodal decomposition 
[9, 10, 18, 19, 20, 21, 25, 26]. Thus domains and topological defects are dynamically deter- 
mined by the nonequilibrium field theory. Domain growth in the Minkowski spacetime ex- 
hibits the Cahn- Allen scaling behavior of classical theory. Further, unstable long wavelength 
modes lose quantum coherence through interaction with stable short wavelength modes and 
show classical correlation, thus achieve classicality. Finally, we apply this nonequilibrium 
field theory to the early expanding universe and study the domain growth and topological 
defects. Domains remain frozen in the infiation era and have scale-invariance. This affects 
the topological density in the early universe. 

The organization of this paper is as follows: In Sec. II equilibrium phase transitions are 
briefly reviewed. In Sec. Ill we discuss the Kibble-Zurek mechanism for the dynamics of 
equilibrium processes. In Sec. IV some quench models are introduced for nonequilibrium 
phase transitions. In Sec.V we introduce nonequilibrium quantum field theory based on the 
functional Schrodinger equation. In Sec. VI the quantum LvN equation is used to solve 
the functional Schrodinger equation. In Sec. VII we find the correlation functions during 
the phase transitions, which determine the domain size. In Sec. VIII nonequilibrium phase 
transitions are studied in an expanding universe. In Sec. IX we show how the classicality of 
classical correlation and decoherence can be achieved for long wavelength modes. 

II. EQUILIBRIUM PHASE TRANSITIONS 

To understand the symmetry breaking or restoration mechanism, we consider a scalar 
field model with the potential 

V{<P) = + (1) 

For a negative the symmetry of the system is broken. However, quantum fiuctuations 
around the true vacuum may restore the symmetry when the temperature is high enough 
so that the thermal energy can overcome the potential energy difference between the true 
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and false vacua. To find the thermal contribution (correction) to the classical potential, 
we divide the field, = 0c + 0/, where 0c is the classical background field and 0/ denotes 
quantum fluctuations around 0c. Then the effective potential is given by 

VM ^ MpV) -l{ml + ^r^)0^ + ^0t (2) 

where p is the density operator, and rriR and Xr are renormalized coupling constants. For 
the broken symmetry — > — m^, one has the potential 

Vei<f>c)-l{T'-TM + ^€ (3) 

The system thus restores the broken symmetry when T > Tc [4]. In other words, the symme- 
try can be spontaneously broken when the temperature drops below the critical temperature 
T 

c- 

In the early stage of evolution after the Big Bang, the Universe would have undergone a 
sequence of phase transitions as the temperature dropped due to the expansion. A possible 
sequence of phase transitions based on particle physics is the GUT phase transition at 
Tc ^ 10^^ — 10^^ Gev, the EW (electroweak) phase transition at Tc ^ 10^ Gev, and Color 
definement/Chiral phase transition at Tc ^ 10^ Mev. Depending on the particle physics 
model, the system produces different type of topological defects [3] . The full symmetry of 
the system is broken to a subgroup after a phase transition. The structure of the vacuum 
manifold Ai, the homotopy group, determines the type of topological defects: domain walls 
for 7To{M.) 7^ oo, strings for 7ii{A4) ^ oo, and monopoles for tt2{M.) ^ oo. 

III. KIBBLE-ZUREK MECHANISM 

Kibble used the principle of causality and the Ginzburg temperature to find the correla- 
tion length of the same phase (domain) in phase transitions [1]. The Ginzburg temperature 
is the one where the thermal energy is comparable to the free energy of a correlated region 

hTc^e{TG)AF{TG), (4) 

so that the field can overcome the potential barrier to jump to other configurations. In 
this case topological defects lose stability. This temperature restricts the size of correlation 
length for stable defects. Topological defects are located along the boundaries of correlated 
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regions. Thus there is one monopole per volume and the density of monopoles is given 
by k/^^, and similarly one string per area and the density by k/^^. 

On the other hand, the Zurek mechanism incorporates the dynamics of equilibrium pro- 
cesses [2]. In the adiabatic cooling (quenching) the equilibrium correlation length and the 
equilibrium relaxation time increase, respectively, as 

e = eo|er", r = ro|e|-'^ (5) 

where e — {Tc — T)/Tc. Here /i and u are model-dependent parameters. Also e is related 
with the quench rate as e = t/rQ. As temperature approaches the critical value (T — >■ Tc or 
t 0), the equilibrium relaxation time becomes sufficiently longer and the process critically 
slows down. However, the correlation length increases indefinitely but the propagation of 
small fluctuations is flnite; sor oc ^/v — > oo. Therefore, there is a time t* when the 
correlation length freezes: \t*\ — T{t*). Prom the above equations the correlation length is 
given by 

m^r^Q^^'^'^. (6) 

Thus the Kibble-Zurek mechanism determines the domain size dynamically in terms of the 
quench rate. 



IV. NONEQUILIBRIUM PHASE TRANSITIONS 

There are many cases to which nonequilibrium phase transitions can be applied but equi- 
librium phase transitions cannot be applied. Matter fields in a rapidly expanding universe 
are such an example, where the Tolman temperature drops as 

m-n^y (7) 

where a is the scale factor of the universe. Another example is the rapid quenching processes 
such as in the rapid cooling of quark-gluon plasma in the Heavy Ion Collision and the liquid 
helium He^ and He^, where domain walls and vortices (strings) may be formed. 

As a field-theoretical model for nonequilibrium phase transitions, we consider the scalar 
field potential [9, 10, 18, 19, 20, 21, 25] 

V{<f>)^lm'm' + ^_^\ (8) 
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where the mass changes signs from m^(— oo) = mf to m^(+oo) = ~"^/- When tq is the 
time scale for the quench, one can distinguish the adiabatic quench process, Am^/rg <^ 1, 
from the rapid quench process, Am^/rg ^ 1. For instance, an analytical model may be 
considered [18] 

= -m^tanhf— ), (9) 

in which for i ^ — oo and — > — for t — > +oo. In the limiting case tq ^ 0, 

one has the instantaneous quench, = -\-w? for i < and = —w? for i > 0. The field 
model can easily be generalized to an expanding universe 

where a is the scale factor of the FRW universe. The Minkowski spacetime is the special 
case of a = 1. In the Hartree-Fock and mean-field approximation, we divide = 0c + 0/ to 
obtain the equations of motion for the classical field 

0e + ^-J>c - V^0e + {m\t) + ^(g + ^ W))</'c = 0, (11) 

and for quantum fluctuations 

h + 3^0/ - V^0/ + (m^(t) + ^<t>l + = 0. (12) 

V. NONEQUILIBRIUM QUANTUM FIELD THEORY 

The functional Schrodingcr-picturc provides the real-time evolution of (time-dependent) 
quantum systems [22, 23]. It is based on the functional Schrodinger equation (in unit of 

i-*(0, t) = Hi<l>, -i—, i)*(0, t), (13) 

where (f) represents a scalar or fermion field. It describes the evolution of wave functionals 
^(0,t) from one spacelike hypersurfacc to another E^. The set of all wave functionals 
constitutes a Hilbert space, which has an inner product on each spacelike hypersurfacc 

The action of operators are defined as 

a(0, tt) |*(0, tt)) ^ 0(0, -i^)*(0, t). (15) 



Now the task is to solve Eq. (13) for the symmetry breaking potential (8). The evolution 
of the wave functional may be found in terms of Green function (kernel or propagator) 

*(x,t) = y^G'(x, t;xo,to)*(xo,to)c?xodto, (16) 

where the wave functional implicitly depends on space through 0. Our stratagem is to 
separate the Hamiltonian into a quadratic part, an exactly solvable one, and a perturbation 
part: 

H{t)^Ho{t) + XHp{t). (17) 

Here Hq includes not only the quadratic potential term but also some contribution from 
non-linear terms a la the Hartree-Fock or mean-field approximation. In this way the approx- 
imation becomes nonperturbative including parts of quantum corrections from the nonlinear 
terms. Then we may introduce a Green function for Sts 

(i^ - ^o(x, t))Go(x, t; x', t') = 5(x - x')5(t - t'), (18) 

and write the wave functional in terms of Go 

^(x, t) = *o(x, t) + xj Go(x, t; X , t')Hp{x', t')^(x , t')di^dt', (19) 

where ^'o is a wave functional for Hq. The wave functional ^ can be put recursively into 
the righthand side of Eq. (19) to result in 

*(1) = *o(l) + A y Go(l, 2)^p(2)*o(2) + A' / / ^0(1, 2)^p(2)Go(2, 3)iyp(3)*o(3) + • • • , 

(20) 

where (i) denotes (xj,fj). 

The general wave functional for Hq can take the Gaussian wave functional 

= 7Vexp[- / (0(x)-0(x,t))( ^ -zS(x,y,t)) 

X (</>(y) - 0(y, ^)) + ^ / ^(x, t) (0(x) - 0(x, t))] . (21) 
It has nonzero expectation values of the field and momentum 

(*o|0(x)|*o) = 0(x,t), (*o|7r(x)|*o) =7f(x,t), (22) 
and also the two-point correlation functions 

(*o|0(x)0(y)|*o) - 0(x,t)0(y,t) + G'(x,y,i), 

(*o|7r(x)7r(y)|*o) = 7f(x, t)7f(y, t) + E(x, y, t). (23) 



VI. CANONICAL METHOD FOR NONEQUILIBRIUM PHASE TRANSITIONS 



Recently, another canonical method, the so-called LvN or invariant method, has been 
developed based on the quantum LvN equation [17, 18, 19, 20, 21] 

^^0{(|>, -i^,t) + [0(0, -i^, t), H{^, -i^, t)] = 0. (24) 

The idea of the LvN method for quantum mechanical systems first exploited by Lewis and 
Riesenfeld [24] is to solve Eq. (24) and find the solution to the Schrodinger equation as 
eigenstates of the operator in Eq. (24). In quantum field theory the wave functional to the 
Schrodinger equation is directly given by the wave functional of the operator 

C>(x, i)*^(x, t) = (^(x)*^(x, t). (25) 

Note that the eigenvalue (/^(x) docs not depend on time, which is a consequence of Eq. (24). 
In particular, for the quadratic Hamiltonian the operator satisfying Eq. (24) can be obtained 
explicitly. This canonical method has an advantage that quantum statistical information 
can be naturally incorporated into the dynamics even for nonequilibrium systems. Another 
merit is that it is a nonperturbative canonical method and can further be applied to fermionic 
and gauge systems. 

We now turn to the potential (8) for nonequilibrium phase transition. As explained in Sec. 
V, we separate the Hamiltonian density Ti into the quadratic part 7io and the perturbation 
part Hp: 

Wp(i) = a='(i#'-i(0'>0'). (26) 
It is convenient to decompose the field into Fourier modes and then to redefine them as 

= ^iMt) + <f>-^m, <f>i~\t) = ^iMt) - '/•-k(i)], (27) 

where (j)^^^ and arc the Fourier-cosine and sine modes, respectively. For simplicity 
reason a compact notation a, (5,- • •, will be used for {k, (±)}. Then the actual Hamiltonian 
of symmetric state modes takes the form 

„ , -i.i^{i\rh^ -I- 
'^2a? " 2 



m - E [h< + ^^imi] + ^ [E + 3 E m] > (28) 



a 



where 

uUt) = m'{t) + ^. (29) 

(Jj 

Thus the Hamiltonian consists of infinite number of coupled anharmonic oscillators, all of 
which depend on time through m^{t). The quadratic part of the Hamiltonian becomes 

Ho{t)=E^s^l + Y^lit)'Pl (30) 



a 



where 

^l(t)-m\t) + ^ + ^j:{<l>l). (31) 
a z ^ 

Note that the Hamiltonian (30) simply consists of decomposed oscillators. 

We first find the Green function Go for Hq and then obtain perturbatively the wave 
functional for the full Hamiltonian. In fact, each mode of the quadratic part Hq can be 
solved exactly in terms the time-dependent creation and annihilation operators satisfying 
the LvN equation. Then the wave functional for Hq is the product of the wave function 
for each mode. The time-dependent creation and annihilation operators of the a mode are 
given by [17, 18, 19, 20, 21] 

We require these operators to satisfy the LvN equation with Ha 

+ [al{t),Ha{t)\ = 0, iQ^aait) + [aa{t),Ha{t)\ = 0. (33) 
Then the auxiliary field ipa satisfies the mean-field equation 

(Pa{t) + 3^0„(t) + nl{t)^a{t) = 0. (34) 

a[t) 

Indeed, these operators satisfy the usual commutation relations at equal times 

[aa(t),aJW] = (35) 
when the Wronskian condition meets 

a^(0>a - <0a) = i- (36) 
In the oscillator representation of {aj^, Sq}, the quadratic part takes the form 

Ut) = y E(<^a + ^>a)«a + '^{^l^a + ^^»a)4«a + {^^^ + , (37) 

^ a. 



and so does the perturbation part 
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EE 

a k=0 



X [ifldl + 2ip*^ippdldp + (fifdf) 



(38) 



The essential point of the LvN method is that the number states of aj^ and Oq, 

Na{t)\na,t)o = dl{t)dait)\na,t)o = na\na,t)o, (39) 

are exact quantum states of the time-dependent Schrodinger equation. The quantum state 
of the field itself is then a product of each mode state. For instance, the Gaussian vacuum 
state of the field is given by 

\0,t)o = l[K,t)o. (40) 



We now find the Green function for Ha 



(41) 



and the Green function for Hq 



G'o(x,t;x',t') = nG'oa(0a,t;0a,O 



(42) 



The wave functional for the full Hamiltonian (28) can be obtained at least perturbatively 
by substituting Eqs. (42) and (38) into Eq. (20). 



VII. DOMAIN GROWTH IN MINKOWSKI SPACETIME 



We now study how the dynamical processes of nonequilibrium phase transitions affect the 
domain growth and topological defect density. The quench models in Sec. IV describe such 
nonequilibrium processes, which can be treated exactly for the free theory and approximately 
for the sclf-intcracting theory. Before the onset of the phase transitions, the mass term 
dominates over the last term from the quantum corrections. It is thus justified to use 
approximately the free theory after the onset of phase transition until it crosses the infiection 
point or spinodal line. This is also true for the quench process lasting for an indefinitely 
long period. We consider first the free theory in the Minkowski spacetime and then discuss 
the effect of nonlinear terms. 
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The free theory for the quench models is provided by the potential (8), where A = and 
m^(i) changes signs either instantaneously or for a finite period. In the Minkowski spacetime 
we can apply the formalism in Sec. VI simply by letting a — 1. Before the phase transition 
(m^ = m^), all the modes are stable and oscillate around the true vacuum: 

1 



+ m?. 



The two-point correlation function is the Green function at equal times 

Go(x, x', t) = (0(x, i)0(x', t))o = Go(x, t; x', t) 



(43) 



(44) 



with respect to the Gaussian vacuum or thermal equilibrium. The two-point thermal corre- 
lation function has the well-known form 



1 miKi{mi\x-y^, 
GioT(y,x,i) = — 1 1 h 



rrii 



47r2 |x-y| 27ry|x-y|2 + m2 

OO 

X Y: i^iK^|x-y|2 + (/?n)2], 



(45) 



n=l 



in terms of the modified Bessel function Ki. However, after the phase transition (m^ — 
—irij), the true vacuum becomes the false vacuum. 

First, in the instantaneous quench model, long wavelength modes {k < rrif) become 
unstable and exponentially grow as 

1 



^/Ok(^) 



I 



sinh(a;y^ki) -I- cosh(a;/kt) , {pf\L — \/"^/ — k^); (46) 



whereas short wavelength modes {k > rrif) are still stable and have the solutions 



<^/ok(t) = 



. <^ik . 



sin(a;/k^) + cos(a;/kt) , (^^/k = ^^k^ - m^j. 



(47) 



The two-point thermal correlation function (r = |x— x'|) is dominated by the long wavelength 
modes and is given by 



G'/oT(r,t) ~G'or(0,t)- 



sml A/^r 



2f ■ ) f irifr- 



exp 



^ 8t 



(48) 



Prom Eq. (48) follows the famous Cahn- Allen scahng relation for the domain size 



V 



(49) 
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This implies that domains grow according to Eq. (49) until most of long wavelength modes 
cross the inflection point and sample over the true vacuum. 

Second, we consider a finite quench model with a finite quench period. Such a model is 
mimicked by the mass term 

m^{t) ^ml-ml tanh(^) . (50) 

The modes with k < irif — ^JruQ — mi become unstable after the phase transition and 
exponentially grow, but those with k > mj are still stable and oscillate around the false 
vacuum. Long after the phase transition, the unstable long wavelength modes have the 
asymptotic solution 

Cn,, - /^k (^2 k2)l/2, (^2,1,2)1/2, 

whereas the stable short wavelength modes oscillate around the false vacuum. Here /Xk and 

z/k depend on the quench process and satisfy the relation |/XkP ~ |i^kP = 1- The correlation 
function is then approximately given by 

G,„(x,,x'.*)^^ (52) 

The coefficients /ik are obtained from the exact solutions of Eq. (34) in terms of the hyperge- 
ometric function [18]. Using the exact solutions, we obtain the two-point thermal correlation 
function during the the quench (— r < i < r) 



sin( ) 
V mo / 



^2 



G„,T(r,i) G^,r{0,t) ^° ^ exp(--^). (53) 



mo mo 

Now the Cahn- Allen scahng relation reads 



^.W = 2e)'- (54, 



• 



After the quench (t ^ r), the two-point thermal correlation function is given by 



sm 



Gfo,j.(r,«) ^ G/„.T(0.f) ""''-V|''^ xp(-!!^), (55) 

where 



T 



3 



t-t--[Ci3)-l](u;l^ + u;i^). (56) 
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Then the Cahn- Allen scaling relation for the domain size takes another form 



\ -■ (57) 



The power of the scaling relation is the same as the instantaneous quench, except for a 
time-lag proportional to the cube of the quench duration r. Surprisingly, the correlation 
function has the pole structure at 

Cbf,\,T^n (n = 1,2,3, •• •). (58) 

This implies larger domains for certain quench rates r [18]. 

Finally, we discuss the effect of nonlinear term to the domain growth. As all wave 
functionals of are readily found, Eq. (20) leads to the wave functional beyond the 
Hartree approximation. Putting the perturbation terms (38) into Eq. (20), we first find the 
wave functional of the form 

oo 

* = ^(o) + ^A"*J"\ (59) 

n=l 

and then calculate the Green function using the wave functional (59). An important result 
is that the correlation length beyond the Hartree approximation has an additional factor 



e(t) = V2^^CD(t). (60) 

This factor from the nth order contribution is a consequence of multiple scattering among 
different unstable modes [21]. The physical implication of the nonlinear effects is that the 
correlation length increases by (2n + 1)^^^, where n is the order of quantum contributions 
which depends the time for crossing the spinodal line, i.e., the period for the field rolling 
from the false vacuum into the true one. 

VIII. TOPOLOGICAL DEFECT DENSITY IN AN EXPANDING UNIVERSE 

As the unverse expands, the temperature drops and thus provides a cooling process in a 
natural way. An adiabatically expanding universe has the Tolman temperature 

m = T,^, (61) 

a{t) 
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where Tc is the critical temperature for the phase transition and Oc is the size at this moment. 
Then the mass term for the second order phase transition may take the form 

We now consider the free theory with the mass term (62) in the expanding universe. The 
Hamiltonian for the free field takes the form 

^0(t)=E^ + ^^aW'^a, (63) 

where 

-^W = -T.^ + ^. (64) 

To find exact quantum states, we follow the formalism in Sec. VI. The time-dependent 
creation and annihilation operators (32) are determined by the auxiliary variables now sat- 
isfying the equation 

(Pa{t) + 3^0a(i) + ^l{t)iPa{t) = 0. (65) 

Rewriting the auxiliary variables as (pa = a~^/'^Vai Eq. (65) can be written in a canonical 
form 

Ut) + (rM\t) + ^l^^)v^{t) = 0, (66) 

where 

In particular, we are interested in the unstable long wavelength modes with 



k < ^^M^a^ - T^al (68) 
These long wavelength modes have both growing and damping solutions of the form 



Clc 



-/®« + -^e-/®«, (69) 



where 

Here ciq, and C2a must satisfy cl^C2a — ciaC^a = 1 to guarantee the Wronskian condition 
(36). Assuming a slow variation of 0a, we obtain the adiabatic (WKB) solution with the 
exponent 

/e„w = /MW-f/^^ + .... (71) 
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Then the adiabatic solutions lead to the correlation function 

1 r d^k 



G/o(x,t;x',t) ~ - / 

(X J 



ik-(x-x') k'lkP 



xexp(2/M-.7J-,) 



(72) 



Compared with the correlation function (52) in the Minkowski spacetime, the correlation 
function (72) has an overall decreasing factor from the expansion of the universe. 

Further, after integrating over momenta, the Cahn- Allen scaling relation for the domain size 
in the comoving frame is given by 

=(«/.! AS?)"- p^' 

where to is the time for most of long wavelength modes to undergo rolling motion from the 
false vacuum to the true vacuum. For instance, during the inflation period with a — e^°*, 
the scahng relation becomes 

1/2 

I I — fi """^^ 1 

-MoHo 

where 



= + ^Hl (75) 

The physical implication of Eq. (74) is that the domain size in the comoving frame is frozen, 
Coit) ~ ^D{t') during the inflation period, in contrast with the Minkowski spacetime where 
domains increase due to spinodal instability. Thus there is no Cahn- Allen scaling relation 
and domains remain scale-invariant in the comoving frame. However, the physical size of 
domains increases in proportion to the scale factor of the expanding universe 

Cpoit) = amD. (76) 

Of course, the higher order contributions beyond the Hartree approximation increase both 
the comoving domain size and physical domain size by additional factor \/2n + 1 for the 
n-th order. 



IX. CLASSICALITY OF QUANTUM PHASE TRANSITIONS 

We investigated phase transitions within the framework of quantum field theory whereas 
phase transitions are studied in classical theory in most literature. Though some of char- 
acteristic features of phase transitions are well described and tested by classical theory, 
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quantum theory is more fundamental than classical theory. So there should be some mech- 
anism for the emergence of classical features from quantum phase transitions. In more 
general context than phase transitions, this is known as the quantum-to-classical transition 
[27]. In the previous section we have shown that long wavelength modes play a key role in 
the dynamical process of phase transitions and thus determine the growth of domains. In 
classical theory long wavelength modes constitute just an order parameter, which describes 
the dynamical process. For an order parameter to achieve classicality, the quantum states 
of long wavelength modes must lose quantum coherence and be classically correlated along 
their classical trajectories [19, 28] . 

To answer this question we consider a simple model motivated by the potential (8). An 
exactly solvable model for a long wavelength mode coupled to a short wavelength mode is 
provided by the quadratic Hamiltonian [19] 

Hit) = + ^u;lit)(t>l + ^4>l + \ujl{t)(t>l + \M2. (77) 

Here 0i and 02 are the long and short wavelength modes, respectively. The former becomes 
unstable but the latter remains stable even after the phase transition. The time-dependent 
Schrodinger equation has a Gaussianal wave function of the form 

^o(0i,02,t) = iV(t)exp[-{Ai(t)0? + \B{t)M2 + ^2(t)02}]- (78) 
Here N is the normalization constant and the coefficients of the exponent are given by 

^^^'^--^2^' ^^^'^--^2^' ^^'^-'umum ' ^^^^ 

where 

-W-HW-^X^Siif)l-W^o. (80) 

-W+NW + A^(i:fgf)>.W=0. (81, 

To model the instantaneous quench in Sec. VII, we assume ujfit) = —uf and a;|(t) = (x;|. 
In the weak coupling limit (A -C Cui,uj2), we find the approximate solutions 

ui{t) = cie^i*, U2{t) = C2e-^^^*, (82) 



where 



02 = ^2 + ^(^^Y + ^(^- - 7^) i^^Y- (83) 
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Then the reduced density matrix for the long wavelength mode takes the form 

pr{<P'i, <Pi) = n*nJ exp[-rc0i,c - rA0i,A - rM0i,c0i,A], 



(84) 



where = (0i + 0i)/2 and 0i,a = (0i — 0i)/2. To order A^, the coefficients of the 
exponent are approximated by 



Ta = 

rjw = 



a;2(a;i + a;|) 



1 + 



A^6<;2 



— 2ia;i 



1 + 



2odi 



+ O2] 
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3oo^ 



(85) 
(86) 
(87) 



2Ql{Q\ + Qlf {Ql + Qiyy^ 2Ql 4odP 
Roughly speaking, a system is classically correlated when its wave functions are peaked 
along classical trajectories or the contours of the Wigner function is close to classical ones. 
And it decoheres when each trajectory loses quantum coherence with its neighbors. Quan- 
tum decoherence is realized when the diagonal term (pc of the density matrix dominates over 
the off-diagonal term 0a- A more precise measure of quantum decoherence and classical cor- 
relation was introduced by Morikawa [29] . According to his measure, quantum decoherence 
is given by 

r- _ 1 pTc _ A 50)2(0)1 — 002) 
"^"^ ~2^T\~2\ 2ooi{odf + odiy 

and classical correlation by 



(88) 



Sec = 



-r2-r2 



5X' 



\UJl — UJ2\ 



(89) 



The system loses quantum coherence (decoheres) when Sqd <C 1 and recovers classical 
correlation when Sec ^ 1- 

We therefore find that the unstable long wavelength mode becomes classical, gaining 
both quantum decoherence and classical correlation due to the coupling to a stable short 
wavelength mode. In field models, the infinitely large number of short wavelength modes 
provides an environment to unstable long wavelength modes and very quickly make them 
classical. Quite recently, field models have been studied for the quadratic, Yukawa, and 
linear couplings to the environment, where quantum decoherence of long wavelength modes 
is observed and the decoherence time is found to be very short compared with other time 
scales [28]. In this way long wavelength modes become classical immediately after phase 
transitions, from which appears a classical order parameter [19]. 
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X. CONCLUSION 



We critically reviewed the nonequilibrium quantum phase transitions by employing the 
canonical method called LvN (Liouville-von Neumann) method. The nonequilibrium phase 
transitions are characterized by the quench time scale shorter than the thermal relaxation 
time scale, so that they cannot have time to relax to thermal equilibrium and evolve out 
of equilibrium. As such nonequilibrium phase transitions we have considered the quenched 
second order phase transitions modelled by the 0^-theory with an explicitly time-dependent 
mass that changes sign during the quench process. This model is a quantum model for the 
classical Ginzburg- Landau theory for spinodal decomposition. 

We used the LvN method to study such time-dependent quantum field theory. The 
wave functionals of the Schrodingcr equation carry all the quantum information of the sys- 
tem evolving from initial data such as thermal equilibrium or Gaussian vacuum to final 
nonequilibrium states. The essential idea of the LvN method is first to solve the LvN equa- 
tion in terms of whose operators exact wave functionals are found for the time-dependent 
Schrodinger equation. In this LvN method it is easy to incorporate the thermal equilibrium 
because density operators are ready to be found. In particular, the exact quantum states 
and the Green function can be found explicitly for a quadratic potential. So our stratagem 
is to separate the Hamiltonian into a quadratic part of the Hartree-Fock type and a pertur- 
bation. In the momentum space and the oscillator representation, the perturbation part of 
our model consists only of quartic terms of the creation and/ or annihilation operators. By 
using the Green function for the quadratic part and solving perturbatively the perturbation 
part, we are able to find higher order contributions to the wave functional. This method 
provides to go beyond the Hartree-Fock approximation. 

As a consequence of the direct application of the LvN approach to the quenched second 
order phase transitions, we obtained some interesting results such as growth of domain sizes, 
topological density and decoherence of order parameter. In the instantaneous quench model, 
the domain sizes scale according to the Gahn- Allen scahng relation which has been observed 
in classical phase transitions. In the finite quench model, the domain sizes have a different 
scaling behavior during the quench. Surprisingly, the correlation function after the quench 
exhibits resonance for certain quench rates. Further, the nonlinear effect leads to a multiple- 
scaling relation beyond the Hartree-Fock approximation. Another interesting phenomenon 
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is the classicality of phase transitions due to instabihty of long wavelength modes. The 
long wavelength modes become unstable, exponentially grow due to spinodal instability and 
thus exhibit classical correlation. The long wavelength modes therefore achieve classicality 
through the coupling to a large number of stable short wavelength modes (environment). It 
is found that the topological defect density can be reduced by factors {2n + lY^"^, which is 
a consequence of multiple scattering at the higher order beyond the Hartree-Fock approxi- 
mation. 

Quantum fields in the early unverse provide a natural framework for nonequilibrium phase 
transitions as the adiabatically expanding universe has the Tolman temperature dropping 
inversely proportion to the scale factor of the universe. It is widely believed that the universe 
would undergo a sequence of phase transitions. However, the phase transitions may not be 
properly treated in finite temperature field theory due to the adiabatic expansion. We 
apply the LvN method to the quenched second order phase transition in an expanding FRW 
universe. The damping effect of field due to the expansion of the universe dramatically 
changes the domain growth. In fact, domains in the comoving frame remain frozen and 
show scale-invariant behavior in the inflation era. This result strongly contrasts with the 
Cahn-AUen scaling relation in the Minkowski spacctimc, according to which domains grow 
as a power of time due to spinodal instability. Thus domains growing through dynamical 
processes suppress formation of topological defects since topological defects form on the 
boundaries of domains. However, the physical sizes of domains increase in proportion to the 
scale factor of the universe. This implies that topological defects may not be suppressed 
through dynamical processes of nonequilibrium phase transitions but through the inflation 
of the universe. 
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